Vortex lines of the electromagnetic field 
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Relativistic definition of the phase of the electromagnetic field, involving two Lorentz invariants, 
based on the Riemann-Silberstein vector is adopted to extend our previous study [I. Bialynicki- 
Birula, Z. Bialynicka-Birula and C. Sliwa, Phys. Rev. A 61, 032110 (2000)] of the motion of vortex 
lines embedded in the solutions of wave equations from Schrodinger wave mechanics to Maxwell 
theory. It is shown that time evolution of vortex lines has universal features; in Maxwell theory 
it is very similar to that in Schrodinger wave mechanics. Connection with some early work on 
geometrodynamics is established. Simple examples of solutions of Maxwell equations with embedded 
vortex lines are given. Vortex lines in Laguerre-Gaussian beams are treated in some detail. 



PACS numbers: 03.50.De, 42.25.-p, 03.65. Vf, 41.20.Jb 



I. INTRODUCTION 



The physical significance of the singularities of the phase of quantum mechanical wave functions has been recognized 
by Dirac in his work on magnetic monopoles [lj . The hydrodynamic formulation of the Schrodinger theory discovered 
by Madelung 2] provided a vivid interpretation of the lines in space where the phase is singular. These are simply 
the vortex lines in the flow of the probability fluid. The velocity field v(r, t) of this fluid, defined in terms of the 
probability current j, is equal to the gradient of the phase S of the wave function ip = Rexp(iS/h), 



CN 1 v = - = tt~- T7T2 = ■ t 1 ) 
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Therefore, the flow is strictly irrotational in the bulk; vorticity may live only on the lines of singularities of the phase. 
Regular wave functions may have a singular phase only where the wave function vanishes, i.e. where 3hp = and 
Ssip = 0. These two equations define two surfaces in space whose intersection determines the position of vortex lines. 
However, the vanishing of the wave function is the necessary but not the sufficient condition for the existence of vortex 
lines. They exist only if the circulation around the line where the wave function vanishes is different from zero. The 



' univaluedness of the wave function requires the quantization of the circulation 

<j> dl-v = 2imh/m. (2) 



The importance of this condition in the hydrodynamic formulation of wave mechanics has been elucidated for the first 
time by Takabayasi [3|. If Eq. © holds for every closed contour, we may recover the phase S (modulo 2itK) from v 
up to a global, constant phase with the help of the formula 

S(r) =m f dl v. (3) 

Early studies of vortex lines were restricted to wave mechanics but Nye and Berry 0, 0, IE have shown that phase 
singularities or wavefront dislocations play an important role not only in wave mechanics but in all wave theories. A 
general review of phase singularities in wave fields has been recently given by Dennis 0, . There is a substantial 
overlap of concepts (but not of the results) between our work and the works of Berry, Nye and Dennis. While they 
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concentrate mostly on the stationary vortex lines that are found in monochromatic fields, we emphasize the time 
evolution. 

More recently, the study of phase singularities and vortices in optics has evolved into a separate area of research, 
both theoretical and experimental, called singular optics. A recent review of this field is given in Ref. [13 • 

In order to find a natural generalization of Eq. ft}, we need a replacement for the wave function ip in electromag- 
netism. A suitable object appears in the complex form of Maxwell equations known already to Riemann \H\ and 
investigated more closely by Silberstein [T^ at the beginning of the last century. In this formulation the electric and 
magnetic field vectors are replaced by a single complex vector F that we proposed to call the Riemann-Silberstein 
(RS) vector [HQ 

F = (E + iB)/V2. (4) 

Maxwell equations in free space written in terms of F read (c = 1) 

id t F = V x F, (5a) 
V-F = 0. (5b) 

The analogy between Eq . Il5;ill and the Schrodinger wave equation is so close that one is lead to treat F as the 
photon wave function [l4| and apply similar methods to analyze the vortex lines and their motion as we have done in 
Refs. [H03 in nonrelativistic wave mechanics. There is, however, an important difference that requires an extension 
of our previous methods: the RS vector has three components instead of one. Thus, there are three independent 
phases tpx, ip2, <f3 — one for each component and it is not clear which combination of these phases should be treated 
as an overall phase of the electromagnetic field. 

In the case of the the Schrodinger wave function, the information about the phase S of the wave function is stored 
in the velocity field v = VS/m. Hence, one may try to find the proper definition of the phase of the electromagnetic 
field by introducing first the counterpart of Eq. and then use the velocity field to reconstruct the phase. The 
natural generalization of the definition Q is (in dimensionless form) 

„ = 1 Ek(F^F k -(VF*)F k ) 

However, as has been noticed already by Takabayasi in his study of the hydrodynamic formulation of wave mechanics 
of spinning particles |17| . this generalization does not work. For a multicomponent field the velocity defined in this 
way cannot be used to reconstruct the phase because, in general, V x v does not vanish. Even though one can still 
give a hydrodynamic interpretation of Maxwell theory based on the formula ft}, the simplicity of the scalar case is 
completely lost |Tsj . 

In the present paper, the phase of the electromagnetic field and the vortex lines associated with this phase are 
defined in terms of the square of the Riemann-Silberstein vector. Since F 2 is a sum of two electromagnetic invariants, 
the structure of phase singularities associated with F 2 is relativistically invariant. This definition of the phase turns 
out to be equivalent (provided F obeys Maxwell equations) to the one used in the classic papers on geometrodynamics 

Despite the fact that F 2 does not obey any simple wave equation, the time evolution of the vortices exhibits all 
the typical features found before by us for the Schrodinger equation. During the time evolution governed by Maxwell 
equations vortex lines are created and annihilated at a point or in pairs and undergo vortex reconnections. 



II. GEOMETRODYNAMICS AND THE PHASE OF THE ELECTROMAGNETIC FIELD 

In nonrelativistic wave mechanics the phase of the wave function can be obtained from its modulus provided 
we also assume that the wave function obeys the Schrodinger equation. As a matter of fact it was shown by E. 
Feenberg [l^ that to determine the phase from the modulus it is sufficient that the wave function obeys some wave 
equation that leads to conservation of the probability, i.e. to continuity equation. A similar reasoning applied to the 
electromagnetic field also enables one to determine the (properly defined) phase of this field. This discovery has been 
made by Rainich pp| in connection with the problem of the reconstruction of the electromagnetic field from purely 
geometric quantities in general relativity. Independently, although much later, this problem was solved by Wheeler 
and coworkers plll22l I23L |24| in the context of geometrodynamics. 

Very briefly, the reconstruction of the electromagnetic field from geometry may be described as follows. The Einstein 
equations 



(7) 
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enable one to determine the energy-momentum tensor T^„ of the electromagnetic field from the Einstein tensor 
fl™ — g^R/2 that is made of the metric tensor and its derivatives. However, the knowledge of the energy-momentum 
tensor alone is not sufficient to determine completely the electromagnetic field. This is best seen from the formulas 
for the components of this tensor expressed in terms of the RS vector: 

Too = F* ■ F, (8a) 
T oi = e ijk F*F k /i, (8b) 
T tj = -FTF j -FfF t +5 ij F*-F. (8c) 

All components of the energy-momentum tensor are invariant under the common change of the phase of all three 
components of the RS vector — the duality transformation 

E' = E cos <p — B sin ip, (9a) 
B' = E sin (f + B cos (p.. (9b) 

Therefore, the overall phase cannot be determined from the energy-momentum tensor. Note, that in contrast to the 
situation in quantum mechanics, even the global, constant phase of F has a direct physical meaning. It controls the 
relative contribution to the energy-momentum tensor from the electric and the magnetic parts. The duality rotations 
@ with a constant value of ip leave the free Maxwell equations unchanged. However, a phase varying in space and/or 
time would modify the Maxwell equations. That is the reason why the Rainich construction works. Namely, he 
has shown that if one assumes that the electromagnetic field obeys Maxwell equations, the phase of the field may 
be extracted from T M „. For this purpose he introduced the following four- vector built from the components of the 
energy-momentum tensor and its derivatives 

w» = — — :z p (io) 

U T 

and used the line integral of w* 1 to reconstruct the phase. 

Our proposal, how to define the phase of the electromagnetic field is much simpler and yet it turns out to be 
completely equivalent to the definition given by Rainich. We shall define the phase of the electromagnetic field tp(x) 
as half of the phase of the square of the RS vector 

F 2 {x) = e 2lv ^\F 2 {x)\. (11) 

In full analogy with Eq. of nonrclativistic wave mechanics, we define a "velocity" four-vector m m as 

_ {F 2 )*df,F 2 - F 2 d tl {F 2 ^* 

Since F 2 is a complex sum of two electromagnetic invariants 



= d^ix). (12) 



F 2 = S + iV = \{E 2 - B 2 ) + iE-B, (13) 



Ufj, is a true relativistic four-vector 



u »- 2(S 2 +V 2 ) ' 1 ' 

This vector has the same denominator (up to a factor of 2 that scales both the numerator and the denominator) as 
the vector u> M defined by Eq. (|1(J|> since TJT T a — 4(S 2 + T' 2 ). However, in general, the numerators of vectors and 
are different. They do become equal when the electromagnetic field obeys the Maxwell equations. The proof is 
straightforward but rather tedious and will not be presented here. 

In our formulation, the square of the RS vector plays the role of the wave function ip. Vortex lines are to be found 
at the intersection of the 5 = and V = surfaces. As in the case of the Schrodinger wave function, at all points 
where F 2 does not vanish, the vector is by construction a pure gradient 

u^(x) = d^ip(x). (15) 

Therefore, one may recover the phase of F by the following line integral 

ip(x) = /"Vt^O- (16) 
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Since the RS vector is univalued, the phases obtained by choosing different paths connecting the points x and Xq 
may differ only by a multiple of 2ix. In other words, the vorticity associated with w M (or with uj„ in the Rainich 
construction) must be quantized 

jd^u^) = 27m. (17) 

The phase defined by Eq. (|16|) is determined up to a global phase ipo: the value of ip(x) at the lower limit xq of the 
integral. The value of tpo cannot be obtained from the energy-momentum tensor. 

Under duality rotations © when ip varies from to 2ir, the vector E' at each spacetime point draws an ellipse 
in the E — B plane. The same ellipse is drawn by the vector B' . These ellipses become circles on each vortex line 
since then the vectors E and B are orthogonal and of equal length. This property lead Berry and Dennis § to name 
the vortex lines associated with the square of a complex vector field the C (circle) lines in their general classification 
scheme of phase singularities. 

The denominator in Eq. (|14f> may be also expressed in the form 

S 2 +V 2 = ^ E2+ 2 B2 ^j -(Ex Bf. (18) 

Therefore, the vanishing of 5 2 + V 2 at a point also means that the electromagnetic field at this point is pure radiation: 
the energy density and the Poynting vector form a null four-vector. One may say that on vortex lines the energy of 
the electromagnetic field moves locally with the speed of light. We would like to emphasize that the velocity of the 
energy flow of the electromagnetic field is not correlated with the vector u^. Even the geometric properties of the 
Poynting vector and the space part of are different. Since S is a scalar and V is a pseudoscalar, the vector 
is a pseudovector. In the simplest case of a constant electromagnetic field the Poynting vector is E x B, while the 
vector Up vanishes identically. There does not seem to exist a physical quantity whose flow can be identified with u^. 
In this respect the situation is quite different from nonrelativistic wave mechanics where the gradient of the phase 
determines the velocity of the probability flow. 



III. SIMPLE EXAMPLES OF VORTEX LINES 



The analogy between the phase of wave function and the phase of the electromagnetic field is not exact. Unlike 
the Schrodinger wave function, the electromagnetic field does not have to vanish identically along the lines where the 
phase is singular. It is only necessary that the field is null i.e. the two invariants S and V vanish. Still, we believe 
that the lines along which the field is null deserve the name of vortex lines. 

The time evolution of the vortex lines embedded in the solutions of the Maxwell equations is quite similar to the 
evolution of such lines embedded in the solutions of the Schrodinger equation. The simplest examples of solutions with 
vortex lines can again be found among the polynomial functions. Such functions may be viewed as long wavelength 
expansions and were found to be very useful in the study of vortex solutions of the Schrodinger equation |l5l IT^ | and 
the Helmholtz equation 0, l25j . Alternatively, these polynomial solutions may be viewed as local approximations to 
the full solution, valid close to the vortex lines under study. In this case one may imagine that in the exact solution 
the polynomial is multiplied by some slowly varying envelope that makes the full solution localized. We shall give at 
the end of this Section an example of such a solution. 

As an illustration of a typical behavior of electromagnetic vortex lines, we present very simple examples of the 
electromagnetic field. The following four fields satisfy the Maxwell equations and possess the vortex structures very 
similar to those found in Schrodinger wave mechanics [TEl Il6| 

F (q) = {y + it,z- a + i(a + t),x + it}, (19a) 

F (b) = {y + t,a-i(z + a + t) lX + it}, (19b) 

F (c) = {2x + y-a + i(z + y + t),z-y + t + i(y-a),-t + it}, (19c) 

= {z 2 + t 2 -iat,a 2 - i(2zt + a 2 + ax),a(y - t)}, (19d) 

where a is a parameter that sets the scale for the vortex configuration. In the first three cases the electromagnetic 
fields are linear functions of the coordinates and in the last case the field is quadratic. In the first case, the two 
invariants are 

S = x 2 + y 2 + (z - a) 2 - a 2 - 2at - 3i 2 , (20a) 
V = 2az + 2t(x + y + z- a) -2a 2 . (20b) 
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The equations 6> = and V = describe a sphere centered at the point (0, 0, a) with the time dependent radius 
\J a 2 + 2at + 3t 2 and a moving plane, respectively. The intersection of these two surfaces is a moving ring shown 
in Fig. n The radius of the sphere decreases for negative values of t until t = —a/3 and then starts increasing. 
The rate of change of the radius exceeds (by a factor of the speed of light showing once again that various 
characteristic features of relativistic fields (like their zeros or maxima) may travel with superluminal speeds without 
violating causality. In this simple example, no change of the topology of vortex line takes place. However, in the three 
remaining cases the topology changes according to the same universal patterns as those found in Schrodinger wave 
mechanics. This universal behaviour of vortex lines is reminiscent of the catastrophe theory |2*7ll2*8j. 

The graphical representation of the motion of the vortex lines in all four cases is straightforward since the equations 
5 = and V = can be solved analytically giving x and y for each value of t as parametric functions of z. In each 
case there are two branches that differ by the sign of the square root. 

x(z,t) = ((a - z)(a + t)± \/a 2 + 2at + 3t 2 ^2t 2 - (a - z) 2 ) /(2t), 

y(z,t) = ((a + t){a-z)T Va 2 + 2at + M 2 y / 2t 2 - (a - z) 2 ) /(2t), (21a) 
x(z,t) = a(a + z)/t — a, 



y(z,t) = -t ± \/t 2 - a 2 \/a 2 - 2at + 2t 2 + 2az - 2zt + z 2 /t, (21b) 
x(z,t) = t+(t 2 ± y/t 4 - 8t 2 (z - a)z - 16(z - a) 2 z 2 ) /(4(z - a)), 

y(z,t) = (a - z)/2+ [t 2 qp y/t 4 - St 2 (z - a)z - I6(z - a) 2 z 2 ^) /(8z), (21c) 



x(z,t) = t±^t 2 - a 2 ^a A + (t 2 + z 2 ) 2 /a 2 , 

y(z,t) = -a- (t 3 + zt(2a + z)) /a 2 . (21d) 

The plots of the functions (b) and (d) show vortex creations and annihilations (Fig.|3and Fig.^J and for the functions 
(c) one obtains vortex reconnections (Fig. [3J). Vortex annihilations occur at t = —a and vortex creations occur at 
t = a. Note that according to the formulas l|21bjl and l|21d|l . at these moments the vortex velocity (dx / dt , dy / dt) 
becomes infinite. 

It is also possible to construct localized, finite energy solutions of Maxwell equation with vortices. We shall give 
just one simple example of such a solution constructed from the following localized solution of the wave equation 

Z(r, t) = {y, x, -t ~ ir}((t + it) 2 - r 2 )" 2 . (22) 

With each vector solution of the wave equation one may associate a solution of Maxwell equations treating the solution 
of the wave equation as a complex counterpart of the Hertz potential. Namely, one may check the RS vector -F(r, t) 
constructed according to the following prescription (2i| 

F(r, t) = V x [id t Z{r, t) + V x Z(r, t)] (23) 

indeed satisfies the Maxwell equations. The square of the vector F has the form 

oo 2(t 2 + 2x 2 + 2y 2 - t 2 ) + z(3x 2 - 3y 2 + 4tr) 

6 ((i + ir) 2 -r 2 )6 • [Zi) 

Since the numerator does not contain the variable z, the vortex lines embedded in this localized solution are straight 
lines parallel to the z axis. Two pairs of such lines are created at t — — r/3 at the points (±2r/3, 0) in the xy plane. 
The four vortex lines move (Fig. |3J until they annihilate in pairs at t = r/3 at the points (0, ±2r/3). The speed of 
each vortex line at the moment of creation and annihilation is infinite, showing very vividly that also for localized 
solutions of Maxwell equations the motion of vortex lines may be superluminal without any limitations. Arbitrarily 
high speed of vortex lines associated with solutions of the relativistic scalar wave equation has already been noted in 

Refs. am. 



IV. VORTEX LINES IN SUPERPOSITIONS OF PLANE WAVES AND IN GAUSSIAN BEAMS 



Solutions of Maxwell equations exhibiting vortex structures may also be obtained with the use of standard building 
blocks — the monochromatic plane waves. A single plane wave is described by a null field since both invariants vanish. 
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FIG. 1: Time evolution of a vortex line in the form of a rotating and expanding ring. All figures in this paper were produced 
with the use of Mathematica |2q| . 




FIG. 2: Time evolution of two vortex lines that at t = a suddenly appear as a straight line and then separate and fly away. 
The same solution for negative times would show a reversed process: the convergence of two vortex line and their annihilation 
at t = —a. 

Therefore, the velocity l|14|) vanishes — a single plane wave has no vortex structure. Also, the sum of two plane waves 
does not have any vortex structure; even though it has a nonvanishing velocity field. However, for three plane waves 
we may have various kinds of vortex structures. As an example, we choose three circularly polarized monochromatic 
waves of the same frequency, handedness, and amplitude, moving in three mutually orthogonal directions. The RS 
vector in this case (up to a constant amplitude) has the form 



F(r, t) = (m + ih)e~ t{t ' l - r) + (n + ii)e-^-^ + (1 + im) e - , ( i -""'' , 



(25) 
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FIG. 3: Time evolution of two vortex lines that at t = are mutually perpendicular and nonintersecting. At the time 
t = a(\/2 — l) 1 / 2 the vortex lines cross and undergo a reconnection. 




FIG. 4: Time evolution of two vortex lines that at t = a suddenly appear along a parabola in the y = a plane (shown as a line 
in the first frame) that first opens up very fast into two symmetric wings that later slowly separate. 

where I, rh, and n are three orthogonal unit vectors, the coordinates are measured in units of the inverse wave vector 
and time is measured in units of inverse angular frequency. The square of this vector vanishes at the points satisfying 
the equation 



e i(l+rh)-r _j_ e i(m+n)-r _|_ e i(n+l)-r _ q 



(26) 



8 




FIG. 5: Time evolution in the xy plane of two pairs of straight vortex lines parallel to the z axis. The evolution is indicated by 
the arrows. Pairs of vortex lines are created at the points denoted by crosses and annihilated at the points denoted by circles. 



■ 




■ 




■ 


15. 






■ 




■ 




■ 


★ 




★ 




★ 


i 






★ 




★ 




★ 




■ 




■ 




"10 




■ 




■ 




■ 






* 




* 




* 




* 




* 




* 




■ 




■ 




■ 








■ 




■ 




■ 












5 1 
















★ 




★ 




★ 


i 


t 




* 




★ 




★ 




■ 




■ 




■ 




■ 




■ 




■ 




-15 


★ 


-10 


★ 


-5 


★ 




★ 


5 


★ 


10 


★ 


15 


■ 




■ 




■ 








■ 




■ 




■ 












-5 
















★ 




★ 




★ 




t 




★ 




★ 




★ 




■ 




■ 




■ 




■ 




■ 




■ 






★ 




★ 




*io 




★ 




★ 




★ 




■ 




■ 




■ 








■ 




■ 




■ 


★ 




★ 




★ 


-15i 






★ 




★ 




★ 




■ 




■ 




■ 




■ 




■ 




■ 





FIG. 6: Positions of vortex lines in the xy plane (measured in units of the inverse wave vector) for the three plane waves moving 
in orthogonal directions. Points marked with different symbols correspond to vortices with opposite circulation. 

It is convenient to chose the coordinate system in such a way that the three basis vectors have the form 

( % \ (fA 

Z = "75 > m = ^2 . n= . (27) 
V 73 / V 7s / \ 73 / 

because then all the vortex lines are parallel to the z axis. The position of the vortex lines in the xy plane is 
determined by Eq. (|26|l . For the choice l|27() of unit vectors this equation has the form (apart from an overall phase- 
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factor exp 

1 + e -i(V3x+y)/V2 + e -i(V3x-y)/V2 _ q ^8) 

The solutions of this equation are 



l~2 2 
x n m = 7r v/ 3 (m + n), y± n = izV2 (±-+m-n), (29) 



where m and n are arbitrary integers. The lattice of vortex lines is shown in Fig. El This example shows that vortex 
lines associated with the phase of the RS vector do not necessarily move; they can also be stationary. 

When one of the polarizations of the three waves, say the last one in Eq. H25[). is opposite, the position of vortex 
lines is determined by a time-dependent equation 

e »(f+fh)-r-2it _|_ e i(rh-h)-r ^ e i(t--h)-r _ q j-^q-j 

In this case it is convenient to choose the orthonormal unit vectors in the form 






™ = ~72 I ' « = I | (31) 
VS J V 

The position of vortex lines in the xy plane is determined by the equation 

I _|_ e -i(S3x+y) I V2+2it + e -i(V3x-y) / ' ^2+2it _ q ^2) 

Thus, in this case the lattice of vortex lines is not stationary but it is moving as a whole with the speed of ^8 /3c in 
the x direction. 

The most interesting case, of course, is a superposition not of a few but of a continuum of plane waves, forming 
a collimated beam. We shall concentrate on the Laguerre-Gaussian beams, in view of their applicability to realistic 
situations (cf., for example [Hi}). We use the representation of these beams in the vector theory as in Refs. pul I32I 
123, IS) but we combine the electric and the magnetic field vectors into the complex RS vector This vector for 
Laguerre-Gaussian beams of circular polarization can be written in the form 

F(x, y, z, t) = e -*M-M (ku, iku, i{d x u + id y u)) . (33) 

The square of the this vector is equal to 

(F(x, y, z, t)) 2 = - e - 2i ^- kz \{d x + id v )u) 2 . (34) 

Note, that the vector F given by Eq. 133JI is not just the analytic signal but the full RS vector as defined by Eq. lf4"| 
whose real part is the electric field and the imaginary part is the magnetic induction. The slowly varying complex 
envelope function u = u(x,y,z) is an arbitrary linear superposition of the functions u nrn (p, cf>, z) defined as (we use 
the notation of Ref. (34|) 



u nm (p, 4>i z ) = Cn 



W(Z) 



exp [—p 2 ] exp 



ikp 2 z 



2(z 2 + z 2 R ) 



(v^)' m W) 



e im< ^g— i(2n+|m|+l) arctan(z/zj?) (35) 



where C nm is the normalization constant, w(z) — wo^/l + (z/ zr) 2 is the z-dependent radius of the beam, p is the 
radial coordinate divided by w(z), L™ is the generalized Laguerre polynomial, and zr = ttwq/X is the Rayleigh length. 
The functions u nm describe the beam with the projection of the orbital angular momentum on the propagation axis 
defined by m. They may be written in the form 



u nm (p, 0, z) = (x ± iy) |m| /„ m (p, z), 

where the upper sign corresponds to the positive values of m. This leads to the following formula 

{d x +id y ){x±iy)\ m \f nm 
= {x± iy)\ m \ ±l p* x {df nm ldp + (1 t 1) Wp). 



(36) 



(37) 
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The velocity (|12t can be obtained by differentiating the phase of the function (d x + id y )u but the expression is quite 
cumbersome. However, it is clear from Eq. (I37[) that the function (d x + id y )u for positive and for negative values of 
m carries m + 1 units of angular momentum in the z direction. Vortex lines defined in terms of the RS vector run 
along the z axis and their vorticity has the strength m + 1. At first, these results seem to be in disagreement with 
the detailed analysis of angular momentum of Laguerre-Gauss beams by Allen, Padget, and Babiker given in Ref. [3^ 
since they have shown that the additional unit of angular momentum is to be added to m or subtracted from m 
depending on the (right or left) polarization of the beam. However, we have broken this symmetry by considering the 
RS vector F and not its complex conjugate. This (arbitrary) choice has fixed the (positive) sign of the polarization. 
With this proviso, our definition of vortex lines in terms of the RS vector leads to the same results as the analysis 
of angular momentum. Each component u nrn has only one vortex line associated with the total angular momentum. 
However, superpositions of several u nm components, depending on their composition, may have additional vortex 
lines. 

The presence of vortex lines in Laguerre-Gaussian beams is due to the definite angular momentum in the direction 
of propagation. The same vortex lines appear also in electromagnetic multipolc fields. In this case the RS vector can 
be written in the form |35| 

F(x,y,z,t) =e-^(k + Vx)jj(kr)(r x V)Y JM (f). (38) 
For the dipole field (J = 1, M = 1) 



„ „. N9 (3 + 4fc 2 r 2 -2fc 4 r 4 + (2fc 2 r 2 -3)cos(2fcr)-6fcrsin(2fcr)) 
(F(x, y, z, t)f = e-^\x + iy) 2 ± V ^ ' - - 11 ■ (39) 

Thus, the dipole field for M = 1 exhibits one vortex line along the z-axis (the direction of the angular momentum 
quantization) with unit vorticity. Higher multipoles will exhibit vortex lines carrying more units of vorticity, depending 
on the value of the z component of the angular momentum. 



V. CONCLUSIONS 



The study presented in this paper fully unifies the description of vortex lines in electromagnetism and in Schrodinger 
wave mechanics. In both cases there is a single complex function of space and time whose phase generically has 
singularities along one-dimensional curves in three-dimensional space — the vortex lines. The velocity four-vector 
associated with the phase of the electromagnetic field plays the same role as the velocity v of the probability fluid 
in wave mechanics. The circulation around each vortex line is quantized in units of 27r. There are two important 
differences. First, the gradient of the electromagnetic phase does not have any obvious dynamical interpretation. 
Second, the electromagnetic field does not vanish identically on vortex lines but only the two relativistic invariants 
vanish and the energy-momentum becomes locally a null four-vector. 

Finally, we would like to mention that in principle one should be able to construct a hydrodynamic form of 
electrodynamics, analogous to the Madelung formulation of wave mechanics. The set of hydrodynamic variables for 
the electromagnetic field would comprise the components of the energy-momentum tensor (only five of them are 
independent, cf., for example |18|) and the velocity vector that carries the information about the phase of the 
RS vector. The quantization condition l|17|) effectively reduces the information contained in to just one scalar 
function giving finally six independent functions. However, we have not found a simple set of equations for these 
hydrodynamic-like variables that would be equivalent to Maxwell theory. 
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